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We study preheating in TV-nation, assuming the Marcenko-Pastur mass distribution, equal energy 
initial conditions at the beginning of inflation and equal axion-matter couplings, where matter is 
taken to be a single, massless bosonic field. By numerical analysis we find that preheating via 
parametric resonance is suppressed, indicating that the old theory of perturbative preheating is 
applicable. While the tensor-to-scalar ratio, the non-Gaussianity parameters and the scalar spectral 
index computed for 7V-flation are similar to those in single field inflation (at least within an ob- 
servationally viable parameter region), our results suggest that the physics of preheating can differ 
significantly from the single field case. 



I. INTRODUCTION 

While inflationary cosmology is becoming a precision 
science, with the advent of recent and upcoming exper- 
iments such as the measurement of the CMBR by the 
Planck satellite [1], the particle physics origin of infla- 
tons still remains unclear. Due to their simplicity, sin- 
gle field models of inflation are considered the most eco- 
nomical explanation of a Gaussian, nearly scale invari- 
ant spectrum of primordial fluctuations, as well as the 
flatness and the large-scale homogeneity of the observed 
universe. Future experiments, however, could change this 
situation and put single field models under pressure. For 
these reasons, there has been a keen interest in build- 
ing multi-field inflationary models, see e.g. the reviews 
@; HI ■ Promising setups of multi-field inflation include 
string- motivated models such as jV-flation [4[, inflation 
from multiple M5-branes [f| and inflation from tachyons 

0- 

Naturally, the large parameter space for couplings, 
masses and initial conditions pertaining to multi-field in- 
flation makes a systematic analysis difficult. An inter- 
esting exception, nevertheless, is TV-flation. W-flation 
is a string motivated implementation of assisted infla- 
tion @, M, H, [HJ EL E3] where a large number of un- 
coupled scalar fields, identified with axions arising from 
KKLT compactification of type IIB string theory, assist 
each other to drive an inflationary phase [H; see also 
One salient feature of this model is 
the possible avoidance of super-Planckian initial values. 
Further, using results from random matrix theory, the 
masses for the axion fields can be shown to conform to 
the Marcenko-Pastur (MP) law [l3j under reasonable ap- 
proximations. The spectrum of the masses is controlled 
by only two parameters, the average mass and a variable 
controlling the ratio between the number of axions and 
the total dimension of the moduli space. This renders Af- 
flation tractable despite the large number of fields. To 
a certain extent, successful inflation is contingent upon 
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the initial conditions, however, the model becomes eas- 
ily tractable by assuming that each field possesses equal 
initial energy. 

Observable cosmological imprints from A/"-flation have 
been computed by several groups. The tensor-to-scalar 
ratio r was calculated by Alabidi and Lyth [l8j and was 
shown to have the same value as in the single field case. 
The non-Gaussianity parameter /jvl wa s computed in 
15, 19] where the deviation from single field models was 
found to be negligible. Unlike r and /nl, the spectral 
index of the curvature perturbation n s depends slightly 
on the model's parameters. [UHjEl] showed that n s is 
smaller (the spectrum being redder) than that found in 
single field models, in agreement with the general discus- 
sion made in [2(3]. These results suggest that the obser- 
vational data predicted by A/"-flation are not drastically 
different from the single-inflaton case. Note however that 
the 5-year WMAP data already excludes some region of 
the parameter space; see plj . 

In this article, we study multi-field preheating, focus- 
ing on A/"-flation as a specific example. To our knowledge, 
a general theory of preheating for multi-field inflationary 
models has not been fully developed. This is in part due 
to the highly non-trivial nature of the string theoretical 
constructions responsible for inflation. However, even at 
the phenomenological level, effects due to multiple infla- 
tons contributing to preheating are largely unexplored in 
contrast to sing le field inflation [22j, |23j, [2J, [25j, [26[ , see 
also [27l. I2H [29l[30l. l3ll. |32[] . For instance, preheating via 
parametric resonance of a matter field might be more effi- 
cient in the presence of multiple inflatons, as indicated in 
Cantor preheating [3L [33l] . Here, a non-periodic variation 
of the matter field's effective mass leads to the dissolution 
of the stability bands and a possible parametric ampli- 
fication of almost all Fourier modes. This expectation 
is based on spectral theory [U [H, [H, H3] , but a quan- 
titative study about the magnitude of the amplification 
with more than two fields is missing [38]. Owing to the 
possible dissolution of the stability bands, one might ex- 
pect that the collective behavior of the fields gives rise to 
efficient particle production after A/"-flation. 

To introduce the notation, let us consider a single field 
model first: the equation of motion for a matter field Xk 
with wavenumber k (assuming the coupling g 2 ip 2 x 2 /2 ) 
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can be written as a Mathieu equation [25j |. 

+ [A - 2q cos(2r)] X k = 0, (1) 

for the comoving matter field — a 3 ^ 2 Xk, where r = 
rricpt is the rescaled time, vn v the inflaton mass, and the 
two resonance parameters are 

g 2 $l k 2 

Here, <J>o is the slowly varying inflaton amplitude, g is 
the inflaton-matter coupling and a the scale factor (see 
also section HV)) . The efficiency of parametric resonance 
is controlled by the resonance parameter q, which needs 
to be large enough (q ^> 1) in order for the resonance ef- 
fect to hold against cosmic expansion. The upper bound 
on the coupling constant g is given by the potential's sta- 
bility condition against quantum gravity effects as well 
as radiative corrections [25|, l39| (unless the potential is 
protected by supersymmetry). This upper bound on g 
and m v ~ 10~ 6 M P (M P is the reduced Planck mass) 
from the COBE normalisation restricts the aforemen- 
tioned resonance parameter q, leaving not much room 
for effective parametric resonance [3{J. One might hope 
to alleviate this fine-tuning in multi-field models. 

For definiteness, we focus on Agnation. We use the 
Marcenko-Pastur mass distribution for the axion masses, 
put forth in [l3| based on random matrix theory, choosing 
the most likely mass distribution, see section [TT1 Further, 
for simplicity, we assume equal energy initial conditions 
at the onset of inflation. Recently, aspects of preheating 
in the context of TV-flation have been considered in [Z(| , 
pointing out the danger of transferring energy preferably 
to hidden sectors instead of standard model particles. 
This reveals an additional need for fine tuning, a possible 
problem for m any string-motivated models of inflation. 
The study in [40j is based entirely on an effective sin- 
gle field description of TV-nation, such as the one above. 
Thus the common lore of parameteric resonance models 
seems to be applicable in this work. Here we take the 
optimistic view that preheating might indeed occur in 
the visible sector. However, we go beyond the single field 
model. We find that at the end of inflation more than 90 
percent of the energy is confined to the lightest 10 per- 
cent of the fields in a very narrow mass range, while the 
remaining heavier fields are already more or less settled 
to the bottom of their potentials. Thus we focus on these 
lightest TV fields during preheating. We also show that 
slow roll is indeed a good approximation, even though 
heavier fields violate the slow roll condition < 1 long 
before preheating starts. The initial values for preheat- 
ing of the crucial N light fields are therefore given by 
their slow roll values at the end of inflation. Since the 
lightest fields are highest up in their potentials, fields 
will join preheating in a staggered fashion. This and the 
still reasonably large number of fields makes a numerical 
treatment mandatory, which we provide. We solve the 



equations of motion for the matter field for various pa- 
rameters following the above prescription and find that 
parametric resonance is less effective for TV-nation. The 
physical reason for the suppression of parametric reso- 
nance is the dephasing of the multiple fields, ft is thus 
clear that effective single field models fail to properly ac- 
count for this effect. We conclude that the old theory of 
perturbative preheating, and not parametric resonance, 
is applicable when many fields of different masses couple 
to a single matter field. 

This article is structured as follows: in section HT1 we 
review TV-flation and its dynamics during slow roll. We 
extend this discussion in section IIII[ where we provide 
an extrapolated slow roll solution that slightly underes- 
timates the inflatons' potential energy during inflation. 
We compare the slow-roll solution with a numerical so- 
lution and argue that even after r\ < f is violated by one 
or more of the heavier axion fields, the overall behavior 
of the inflatons is still well approximated by the slow-roll 
regime, up until preheating commences. In section IfVI 
we set the stage for preheating corresponding to the end 
of slow roll for the effective single field model. To begin, 
we set all axion masses equal to each other and discuss 
the physics of preheating in this particular system; then 
we return to TV-flation, where the Marcenko-Pastur mass 
spectrum is used. Finally, in section [V] we conclude with 
comments and prospects for studying further issues on 
multi-field preheating. In Appendix [X] we give a semi- 
analytic solution that provides an upper bound for the 
inflaton potential, which further supports the observation 
made in section Hill 



II. 7V-FLATION and slow roll 

The action for TV minimally coupled scalar fields, re- 
sponsible for driving an inflationary phase, can be written 
as (see for a review on multi-field inflation) 




( 3 ) 

where we assume canonical kinetic terms. The unper- 
turbed volume expansion rate from an initial hypersur- 
face at t* to a final hypersurface at t c (below we use * 
and c to denote values evaluated at and t c ) is given by 

N(t c ,Q = J Hdt, (4) 

where N is the number of e-folds, H is the Hubble pa- 
rameter and t is cosmic time. 

In TV-nation Q , the TV scalar fields that drive inflation 
are associated with axion fields. All cross-couplings van- 
ish when the periodic potentials are expanded around 
their minima [13j |. Therefore, in the proximity of their 
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minima the fields have a potential of the form 

Af 

W(iputp2,...,(pAf) = ^2Vi(<fi) 



i=l 



El- 



,2,„2 



(5) 



where the fields have been arranged according to the 
magnitude of their masses, namely wn > rtij if i > j. 
TV-flation is a specific realization 1 of assisted inflation 
@, S EH , where the Af scalar fields assist each other in 
driving an inflationary phase. In this manner, individual 
fields do not need to traverse a super-Planckian stretch 
in field space. The spectrum of masses in (|5|), which were 
assumed to be equal in [J] , can be evaluated by means of 
random matrix theory within a context of KKLT moduli 
stabilisation (44|, and was found by Easther and McAl- 
lister to conform to the Marcenko-Pastur (MP) law [l3l |. 
This results in a probability for a given square mass of 



p(m 2 ) 



1 



2nm 2 f3m* 



(m 



2 

max 



m 2 )(m 2 



,) .(6) 



where (3 and m 2 completely describe the distribution. 
Here, m 2 is the average mass squared and {3 controls 
the width and shape of the spectrum. The smallest and 
largest masses are given by 



""max 



= m(l- y/0), 

= m(l + V/3)- 



(7) 

(8) 

) into 



In this paper we split the mass range (m m j„, m mt , 
Af bins, 



(mo, mi), (mi, m 2 ), • • • , (msf-i,rhjs) (9) 
where mo = m rn i n , fhj\r = m max and mj_i < m,, so that 

~ 2 

1 



p[m )dm 



Af' 



1,2, •••,7V". (10) 



We then represent each bin (mj_i,mj) by an inflaton of 
mass m.j. In practice we simply set 

m 2 = (m 2 _ 1 +m 2 )/2, i = 1, 2, • • • ,Af, (11) 

in the numerical computations. Apart from the Af in- 
flatons tp±, ip2, ■ ■ ■ , <pn with masses mi,m2, • • • , mj^, we 
introduce a fiducial inflaton <po with mass mo = m m i„ for 
computational convenience (we shall use this as a clock). 
In [13j, /3 is identified with the number of axions con- 
tributing to inflation divided by the total dimension of 



1 For a different realization of assisted inflation based on M-theory 
from multiple M5-branes see 0] [4l| . See also [42j for another 
approach to random potentials in the landscape. 



the moduli space (Kahler, complex structure and dilaton) 
in a given KKLT compactification of type IIB string the- 
ory. Due to constraints arising from the renormalization 
of Newton's constant Q ~ 1/2 is preferred. Hence, 
we will work with (3 = 1/2 in the following. Further, 
the magnitude of m is constrained by the COBE nor- 
malization (l3l . [45l | , so that there is not much freedom in 
A/"-flation to tune parameters. 

At this point, we introduce a convenient dimensionless 
mass parameter 



X j 



as well as the suitable short-hand notation 



i + Vg 



(12) 



(13) 



A properly normalised probability distribution for the 



variable x = m . /m min is p(x) 



% p(m ); hence ex- 



pectation values with respect to the MP-distribution can 
be evaluated via 



</(*)) 



1 " 

i — l 



P(x)f(x) dx 



27T/3 



y/(£-x){x-i)tt&-dx. (14) 



In section IIIIDI and in appendix [A] we make use of the 
additional notation 



</(*)> 



_ (1 - VP? 

2it(3 



^(t- x ){ x -l)M- dx ,(15) 



where a and b are the limits of integration. When f(x) 
is a polynomial in x, (p~4|) reduces to a hypergeometric 
integral. In particular, 



-1/2 



(1) 



1, 



(1 - VP) 2 



(16) 



Figfflshows a plot of the probability distribution function 
p(x), for 13 = 1/2. 

Initially, we restrict ourselves to the slow roll approxi- 
mation. As we have already seen, Af fields contribute to 
the energy density of the universe through a separable 
potential. In this regime, the dynamics of A/"-flation is as 
follows: first note that the field equations and Friedmann 
equations can be written as 



3# 2 



M = -V! 



W . 



(17) 
(18) 



Here and throughout most of our analysis, we set the 
reduced Planck mass to Mp = (SnG)" 1 ^ 2 = 1. The slow 
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FIG. 1: A plot of the Marcenko-Pastur mass distribution over 
the dimensionless mass variable x = m 2 /m^, in , for = 0.5. 
Note that the MP distribution peaks in the small-mass region 
so that the majority of fields resides there. We use mass 
parameters Xi — m^/rn^^ corresponding to the mass-bins 
chosen according to (J9]), (fTtJl) . (fTTj) . so that xq = 1 < x\ < 
X2 < ■ ■ ■ < xj\f-i < xu < ^ 34. 



roll approximation is valid if the parameters 



£l ~ 2 W 2 ' 
are small (e* <C 1, rji -C 1) and 



V" 
W 



e = £i < 1 

i=l 

holds. The number of e-folds of inflation becomes 

i=l J * i 

and the field equations can be integrated to yield 



(19) 



(20) 



(21) 
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III. EVOLUTION OF INFLATONS AND AN 
EFFECTIVE SINGLE FIELD MODEL DURING 
INFLATION 

In this section we investigate the evolution of the ax- 
ion fields after the slow roll condition is violated for one 
or more of them. In order to study the system's collec- 
tive behavior, it is useful to use an effective single field 
description 0. Given equal energy initial conditions for 
the fields, the slow roll parameter r/j\f of the heaviest field 
will be the first one to become of order unity (46|. Hence, 
prior to this moment we can safely implement an effec- 
tive model composed of a single field a which evolves 
according to an effective potential W e ff(a). After ij^ 
became of order one, the corresponding field ip_\f cannot 
be described by the slow-roll solution. Below we argue 
that in our particular model of TV-flation the whole sys- 
tem is nevertheless well approximated by the slow-roll 
solution, up until the slow-roll parameter e of (f3"2"|) be- 
comes of order one 2 ; during this stage, the contribution 
of the heavy axions is negligible compared to that of the 
lighter fields. This is due to three characteristic features 
of the model: (1) the majority of the axions is distributed 
around the lightest mass in the MP law; (2) the small 
value of the heaviest fields is prescribed by the equal en- 
ergy initial condition; hence, heavy fields provide a small 
contribution from the onset; (3) when the slow-roll con- 
dition r)i < 1 is violated for the heaviest fields the Hubble 
parameter is still very large, resulting in an over-damped 
evolution of the heavy fields. 

Naturally, it is possible to continue using the slow roll 
approximation when the next heaviest field violates slow 
roll and so on and so forth. This regime ends when slow 
roll fails for the effective single field a, after which light 
fields will actually start to evolve faster and preheating 
starts. It is important to note that we can trust our ap- 
proximation up until preheating starts, where possible 
particle production due to non-linear parametric reso- 
nance is our main concern. 



A. Effective single field model 



Notice that this relationship between fields does not cor- 
respond to an attractor solution and predictions of TV- 
flation can depend on the initial conditions. Recall that 
we are assuming equal energy initial conditions, namely 



V* — V*, which can be rewritten as 



rriifj)* = in j (j) 



3^3 



(23) 



A subtle feature of TV-nation is that if the mass spectrum 
is broad (£ 3> 1, corresponding to (3 ~ 1), the heavier 
fields will acquire r/i 1, even as inflation continues. 
This is the case even for the preferred value of j3 ~ 1/2, 
corresponding to £ ~ 34. In the next section we show 
how to deal with TV-flation after rji ~ 1 for the heaviest 
inflaton. 



Here we derive the effective single-field model based 
on the slow roll approximation. This provides a lower 
bound to the evolution of the total potential energy. We 
identify the effective inflaton field a as the path-length of 
the trajectory in the TV dimensional field space; namely, 
for TV scalar fields <pi we have Q 



t AT 

E 

* i=l 



ctiifiidt , 



(24) 



Or the slow roll parameter of the effective degree of freedom in 
l32t . 
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with 



<P 2 



(25) 



where the ipi can be computed given the dynamical re- 
lations in (|22|) , which are valid during slow roll, as well 
as the initial conditions in (f2U|) . Note that a — at the 
initial time i*. 
Using 



^0 



(26) 



(where y parameterizes how far the field tpo rolls down 
its potential) and Xi is defined in (fT2"|) , as well as the 
equal energy initial conditions (|23[) . we can rewrite the 
dynamical relations as 



2 y 



(27) 



We are using the fiducial inflaton tpo for computational 
convenience; (po is not one of the Af inflatons driv- 
ing AA-flation (note that there is a vanishing proba- 
bility for to = too = m m i n according to the MP 
law). From the Klein-Gordon equations during slow 
roll along with the Friedmann equation we obtain ipf — 
mfai<pfy Xi /(3W), as well as dy = -(2m%y/V3W)dt, 
with W = (1/2)t77^S 2 E ' 



A' 



1 y Xi . Using these relations 
in (f2"4"|) . we obtain an effective single-field solution (for 
which the subscript / is used), 




1/2 



ds 
s 



-tfo I \J (xs x ) — , 

}y S 



(28) 



where in the last step we used the definition of the MP- 
expectation values from (|14p . Similarly, the correspond- 
ing potential in ([5]) can be computed as 



Wi{y) 



1 

2 

Af 



AT 

„2 *2 a 

"Wo l^V 

i=l 



"w5 2 (y x ) 



(29) 



Equations ([28]) and (f29|) provide an implicit means of 
computing Wi(aj). The number of e- folds can also be 
computed with this approximation. From (|2ip we find, 



Ni(y) 



1 



,„*2 



Af j x . 



A/'' 



(30) 



Let us look at the solution more quantitatively. Since 
we would like to avoid super-Planckian initial conditions 



we assume the lightest possible inflaton to set off from 
the reduced Planck scale, ip^ = 1. Then it follows from 
the equal energy initial conditions that all the other fields 
evolve safely within sub-Planckian scale. Here and in the 
following we use the preferred (3 = 1/2 so that the ratio of 
the heaviest to the lightest mass squared in (fT5|) is about 
£ ~ 34. The e-folding number (|30f> depends linearly on 
the number of inflatons Af. If we take Af = 1500 we get 
N max = Nj(0) « 64.3, which is large enough to solve the 
standard cosmological problems. Note, however, that one 
cannot trust (f2"8")l and (j2"9")l down to y = 0, since slow roll 
ends earlier. Strictly speaking, our effective single field 
solution (with subscripts '/') is only valid as long as the 
slow roll conditions are satisfied, that is until rjjj = 1. 



Using 



and to jv 



this can be written as 



(y x ) 



2£ 

Af^l 2 



(31) 



from which the value of y at rj^f = 1 is found numerically 
as yx « 0.488 for Af = 1500. The number of e-folds at 
this instant is Ni(y^f) ss 55.6 and we see that there is 
still a breadth of inflation to come. If we ignore this fact, 
we could extrapolate 07 up until this effective degree of 
freedom leaves its own slow rolling regime when 



_ 1 (W'j 

Wi 



= 1 . 



This equation can be rewritten as 

2{xy x )=Af^ 2 (y x f , 



(32) 



(33) 



where we used 



1 da! 



A/" 



= E* 



difi 



i=l 

E 

i=l 



«woV 'N ( x v x ) 



(34) 



Equation (|33[) can be numerically solved to obtain y en d ~ 
0.0836 for Af = 1500 so that oi(y end ) m 17.6 and 
Ni{yend) ~ 63.8. At this point inflation comes to an 
end and preheating is about to commence. The poten- 
tial at this instant is Wi{y en d) ~ 0.128777 2 w 1.49mg (see 
Table H. 



B. Implications for non-Gaussianities 

In contrast to single-field inflationary models, in which 
non-Gaussianities (NG) are known to be suppressed, in 
multi-field models there is a possibility that NG may be- 
come large due to the existence of isocurvature pertur- 
bations [43| ; if there is a sudden turn of the trajectory 
in field space as in the case of the curvaton scenario, a 
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-0.2 L 



for M = 
cpt with 



1500 



FIG. 2: Evolution of the axion fields 
(numerical result). The figure shows {pi witn i 
1, 300, 600, 900, 1200, 1500 from the top. The time r = m t is 
in units measured by the fiducial mass mo = m m in- In this 
unit, y — y^f is at r = 7.14 and y = y en d is at r = 11.6. 



conversion of isocurvature perturbations into adiabatic 
ones takes place, giving rise to larger NG. In TV-nation 
the non-linearity parameters characterizing the bi- and 
tri-spectrum were investigated in the horizon crossing ap- 
proximations in [l9[ and it was found that AA-flation is 
indistinguishable from single field inflationary models in 
this limit. Also, incorporating the evolution of pertur- 
bations after horizon crossing, but still within slow roll, 
revealed that additional contributions remain negligible 
[ll| . Hence, NG are expected to be heavily suppressed as 
long as slow roll is considered [4(J. In the present paper 
we investigate the evolution in A/"-flation after the slow 
roll condition is violated for one or more of the axion 
fields and find (see the next sections) that the extrapo- 
lated slow-roll solution remains a good approximation up 
until preheating commences. Thus, in this intermediate 
regime (after slow roll inflation but before preheating), 
NG should also be suppressed; additional NG would be 
due to the evolution of the adiabatic mode after horizon 
crossing; for this to occur, isocurvature modes have to 
source the adiabatic one, but in A/"-flation, the trajectory 
in field space is smooth; therefore, NG should be heavily 
suppressed up until slow roll fails for the effective single 
degree of freedom, i.e. when preheating begins. Dur- 
ing preheating NG may still appear; this requires further 
study, and we hope to come back to this issue in the near 
future. 



ure clearly shows that heavy axions, even (^300, roll down 
the potential rapidly and their oscillation amplitudes are 
much smaller than that of the lightest field <p\. Indeed, 
fields are usually over-damped up until preheating starts. 
Naturally, the lighter fields are expected to be responsi- 
ble for preheating (unless the coupling between the heavy 
fields and a matter field is extremely strong). 

In Table Q] we summarise the values of W and N ob- 
tained by the extrapolated slow-roll solutions (7) and nu- 
merical results, at a = <j(yj^) and o(y eri d)- A compari- 
son between the analytic and numerical computation for 
N = 100,200,400,1500 reveals an agreement, roughly 
within 15%, indicating that the extrapolated slow-roll 
solution (7) is a good approximation up until y e nd- 

In 

Appendix [X] we provide a semi-analytic computation of 
an upper bound to W (solution 77), which further sup- 
ports this observation. How many inflatons still satisfy 
the slow-roll condition rji < 1 at y e nd,7 This can be found 
by comparing the values of W(y en d) / m,Q and the mass pa- 
rameter Xi, since rji = mf/W = Xirn^/W. Numerically, 
we find that 13, 18, 25, 56 lightest fields are still in the 
slow-roll regime at y = y end , for N = 100, 200,400, 1500. 
Analytically, solution 7 yields somewhat smaller values: 
10, 13, 18, 41, respectively. 









solution I 


Numerical 


N 






Wi(y*) 




W{yu) 






Vend 


<Tl(Vend) 


Wi(y end ) 


Ni{y end ) 


W(y end ) 


N(y end ) 


100 


0.964 


0.541 


34.0 


0.758 


34.0 


1.10 




0.502 


3.27 


2.43 


3.67 


2.83 


4.34 


200 


0.879 


2.00 


34.0 


3.72 


34.3 


4.14 




0.331 


5.40 


1.97 


8.00 


2.27 


8.83 


400 


0.762 


4.37 


34.0 


10.9 


34.4 


11.4 




0.211 


8.33 


1.73 


16.6 


1.95 


17.6 


1500 


0.488 


12.86 


34.0 


55.6 


34.4 


56.3 




0.0836 


17.61 


1.49 


63.8 


1.62 


65.1 



TABLE I: Comparison of analytic and numerical solutions 
for the effective single-field values W and N at tj(yu) and 
""(j/end), for the number of inflatons M — 100,200,400, and 
1500. The values of 07 are found using ([28} and the cor- 
responding analytic and numerical values for W/rriQ and N 
are shown. Apart from the conspicuous disagreement in the 
e-folding number TV for small TV, the extrapolated slow-roll 
solutions (I) are relatively in good agreement with the nu- 
merical solutions, roughly within 15% difference. Typically, 
the results of solution (I) slightly underestimate the potential 
W. 



C. Numerical solutions 



Fig[2] shows the time evolution of Af — 1500 axions in 
our setup, namely the MP mass distribution with (3 = 0.5 
and the equal energy initial conditions, obtained numer- 
ically. The plot shows ipi with i = 1,300,600,900,1200 
and 1500. For the initial conditions we used (|2"3")) with 
ifQ = 1 and ip* = 0. Due to these initial values, lighter 
axions evolve from larger values, closer to 1. The fig- 



D. Light axion dominance 

From Fig[2] we can infer that heavy fields lose energy 
quite rapidly and that the later stage of AA-flation is 
driven solely by the light fields. To see this quantita- 
tively let us introduce the ratio of the potential energy 
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larger. This in turn justifies our usage of the approxima- 
tion given by solution /, which corresponds to using slow 
roll for all fields, even if the slow roll condition r\i < 1 is 
violated for the majority of fields; all in all, heavy fields 
do not contribute much to the total potential energy, and 
henceforth inflation. 



Af 


solution I 


Numerical 


100 


62.9 


63.3 


200 


78.4 


80.0 


400 


88.0 


89.9 


1500 


96.5 


97.5 



TABLE II: The ratio of the potential energy carried by the 
lightest axions (10%) with respect to the total potential en- 
ergy, evaluated at the end of slow-roll regime a — a(y en d). 
The table shows both semi-analytic and numerical results. 



FIG. 3: The ratio of potential energies of the light fields to 
all axions, is plotted against A = l/Af and y. The figure 
does not depend on Af. Note that in calculating R{ we are 
using the slow-roll approximation which is not reliable after 
Vend- For larger Af ', the value of y e „d tends to be smaller, and 
hence the domination of light axions tends to be stronger. 



of £ lightest fields to that of all Af fields, defined by 
Wug h t _ ELi Vi 



Hi 



W total J^Zi Vi 



(35) 



Using the extrapolated slow-roll solution (/), this ratio 
can be evaluated as 



Ef =1 (r 



(y x ) 



(36) 



The ratio of the number of light fields to all fields is 
similarly 



I x f 

II * 



(37) 



Figj3] shows Rj(y) versus A and y. Note that this plot 
does not depend on Af (however y^j- and y en d do depend 
on Af). For small y (i.e. at late time), R\ approaches 1 
for even small values of A, indicating that the potential 
energy is dominated by light axions. 

Table [TT] summarizes the amount of the total poten- 
tial energy carried by the lightest 10% of the axions at 
o = cr(yend), for Af = 100,200,400,1500. These values 
are calculated using both, solution I and numerical com- 
putations. The results clearly show that for large Af, 
the potential energy is carried by only a small portion 
of lightest axion fields. Using (|36p. it is not difficult to 
check that this tendency becomes stronger when Af is 



IV. PREHEATING 

Now let us discuss the physics of preheating. In order 
to solve problems of the standard big-bang cosmology 
the number of e-foldings must be large enough; with this 
in mind, we assume TV = 1500 in the following. If we 
forget about its string theoretical origin (which we do 
in this paper), Af can be even larger, giving rise to a 
(harmless) larger number of e-folds. Our choice of Af is 
for definitencss, and also for numerical tractability. 



A. The model of preheating 

We first recall that the axion mass scale is observa- 
tionally constrained. The power spectrum for multi-field 
models is computed by the SN- formalism [48| and for 
A/"-flation it becomes fijl [491 



1 1 

96tt 2 M|, 



E 



96tt 2 M]> 



(38) 



where in the second line the MP distribution and equal 
energy initial conditions have been used. To be specific, if 
we evaluate Vn at t = (corresponding to large scales) 
we obtain 



_Afmfof 
U 96ir 2 M 6 p ^ 



(39) 



Comparing this with V-r ~ 2.3 x 10 9 from WMAP mea- 
surements [49|, [5(1 HH , we find the mass of the lightest 
field 



too ~ 2.4 x 10" 6 M P . 



(40) 
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for TV = 1500, (3 = 0.5 and ip^ = M P . Consequently, the 
average mass is m = m /(l — \/jT) w 8.1 x 10~ 6 A/p . In 
what follows, we assume these values for the mass param- 
eters; all other masses follow via the MP-distribution. 

In the previous section we have seen that the late stage 
of TV-nation is mainly driven by a few light axions; we 
naturally expect that preheating is triggered by these 
lightest axions, coupled to a matter held, at least before 
back reaction becomes important. Thus, in this section 
we focus on the relevant TV = 150 light axions, which 
carry more than 95 % of the total potential energy at 
V = Vend- It is important to note at this point that the 
relevant mass scale for preheating is set by the mass of 
the light fields and not the average mass m. For the mat- 
ter into which the inflatons decay, we consider a massless 
bosonic field x coupled to the axions via the coupling 
\g 2 <PiX 2 > where we assume for simplicity an identical 
coupling constant to each ip^. The model we consider 
is then described by the following Lagrangian, 

i=l ' 



(41) 



The equations of motion are 

ipi + 3Hipi + (mf + g 2 ( X 2 )) <Pi = 0, (42) 

Xk + 3H Xk + ^+S 2 X>i^Xk = 0, (43) 

i i 

+ \(X 2 ) + \9 2 (X 2 )Y,VI ( 44 ) 

i 

where Xk is the mode operator of the matter field and 
(•) is the mode sum over k. We consider the axions and 
gravity as the background, and ignore backreaction from 
the matter field Xk', consequently, (x 2 ) an d (x 2 ) are set 
to zero. 



B. Parametric resonance in the equal-mass case 

Before addressing the more involved preheating sce- 
nario of TV-nation, we discuss a model with TV = 150 
inflatons having the same mass, m, = m. 

In this case, the equal energy initial conditions set the 
same initial values for all TV axions, so that the evolution 
of the TV axions is identical. Neglecting backreaction of 
the matter field we can write the equations of motion as, 



(pi + 3H(pi +m <fi = 0, 
fk 2 

Xk + 3H Xk + (— +.A/W? 



a" 



Xk = 0. 



3H 2 



TV 



($ + m 2 <p1) 



(45) 



Defining ip = v Af(pi , the equations of motion reduce 
to those of the well-understood single field model, yield- 
ing non-perturbative preheating [3U231 H3, HE H2| ■ The 
Klein-Gordon equation for tp reads 



tp + 3Hip = —to tp , 



(46) 



whose solution is approximated during the preheating era 
by 



(p(t) = <j>(t) sin(mi) , 



(47) 



where $(t) = y/8/(y/3mt) [H[ is a slowly decaying ampli- 
tude due to Hubble friction. The corresponding equation 
for a Fourier mode of the matter field reads 



Xk 



' k 2 



2 2 

9 ¥> 



Xk = 0, 



(48) 



where p = k/a is the physical momentum. Due to the 
oscillations of the inflaton field, the mass of the matter 
field becomes time dependent and resonances can occur. 
To see this, introduce q — g 2 ^ 2 /Am 2 , r = mi, — 
2q + k 2 /m 2 a 2 and Xk = a 3 ^ 2 Xk so that |48|) becomes 



d 2 X k 

dT 2 



+ (A k -2qcos{2T))X k = 0, 



(49) 



where we also neglected the term proportional to the 
pressure, — (3/4)(H 2 + 2a/ a). If we ignore the time de- 
pendence of the amplitude $ in q and of A k , Eq. |49|) 
is the Mathieu equation. It is known that parametric 
resonance occurs for wavenumbers k within resonance 
bands (see [25|, [H| for the stability/instability chart). 
This means if k is within the n'th resonance band, the 
corresponding mode increases exponentially 



X k oc 



(50) 



where /^™' ) > is the Floquet index [53j . Physical param- 
eters correspond to the region A k > 2q and, in particular, 
the zero mode k = evolves along the A k = 2q line from 
large q to q ~ 0, as the inflaton amplitude $ decays 
slowly. As it evolves, the system crosses resonance bands 
where exponential particle production takes place. Par- 
ticle production is efficient in the large q (q 3> 1) region, 
broad resonance (or stochastic resonance when expansion 
effects are included). For small q the resonance effect is 
limited as it is not strong enough to hold against red- 
shifting Xk oc a -3 ^ 2 . Then the main concern is whether 
it is possible to have a large enough q in a given model. 
A stringent constraint comes from radiative corrections, 
restricting the value of the coupling to g < 10~ 3 [25l l39l| . 

From the above discussion it can be deduced that 
having many inflatons does not increase q, and reso- 
nance effects are not enhanced. Given that the equa- 
tions of motion lead to the one of the single field model, 

(p(= \/rV ipi) starts oscillating from the single field value 
~ 0.2Afp. Each inflaton ipt oscillates with smaller ampli- 
tude $/vCv", while q is unaltered. In the next subsection, 
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10 20 30 40 50 

FIG. 4: The evolution of the oscillating term ip 2 that drives 
parametric resonance. The horizontal axis is the dimension- 
less time r = mot. Since ip — VTv' fi and tpi are chosen to 
start from Mp — 1, ip 2 starts from M = 150 at r = 0. 

we compare the equal mass case to a broader mass distri- 
bution (MP mass distribution). To this aim, we provide 
numerical plots. In FigHl we show the evolution of the 
oscillating term ip 2 in the equal mass case, where we have 
chosen J\f = 150, m = rn Q = 2.4 x 10~ 6 M P . The initial 
values of tpi are all taken to be Mp = 1, and the initial 
velocities are tpi = 0. In Figf5]the evolution of the mat- 
ter field mode function \k and the comoving occupation 
number of particles, defined by 




are shown for g = 10 3 and k/ai n it = 6.0 x m (corre- 
sponding to a fastest growing mode, see [25|]). Here, 

which reduces to \Jk 2 /a 2 + g 2 f 2 in the present case. The 
condition on the resonance parameter q = g 2 ^ 2 /4m 2 > 
0(1) corresponds to \ip\ 2 > 10~ 5 for g = 10~ 3 , which is 
roughly r < 4000. We are ignoring both, backreaction 
and rescattering effects. The corresponding plots for the 
Marcenko-Pastur distribution are shown in Figs. [5] and 
HO] below. 



C. Numerical results for the Marcenko-Pastur case 

Let us turn to the question of multiple fields whose 
masses obey the MP law. The equations of motion of 
our system are (|42[l. (j43|) . (|44|) ; we are assuming that all 
inflaton fields are coupled to the same matter field with 
identical strength g 2 , and we consider only M = 150 ax- 
ions since the heavier 90% of the axions are negligible in 
the later stage of jV-flation. We also ignore backreaction, 
meaning we set (x 2 ) = (x 2 ) = (which, in the end, is 
justified since amplification of the matter field is found 




FIG. 5: The evolution of (the real part of) the mode Xfe 
(a) and the occupation number nt (b). The initial con- 
ditions for \k ar e set by the positive-frequency solution at 
r = 13.5, when the slow-roll conditions break down. The 
coupling constant is g = 10 -3 and the wavenumber is cho- 
sen as k/am = 6.0 at r = 13.5. One can see amplification 
due to typical stochastic resonance, i.e. the overall amplitude 
grows exponentially while there are occasional decreases of 
the amplitude. We are ignoring backreaction so that reso- 
nances are present until q ~ 0(1), corresponding to r ~ 4000 
for g = 10 -3 . Backreaction from the matter field shuts off 
resonances earlier. 

to be suppressed). The initial values and the velocities 
of the axions at the onset of preheating y — y en d corre- 
sponds to the extrapolated slow-roll solution, discussed 
m Section [1*1*11 



Vi(v«ud = <Po\ ^ , (53) 
y %i 

The initial conditions for the matter field is set by the 
positive frequency mode funtion, Xk(t) = a 3 / 2 Xk(t) — 
e - i ^ fc (t-r c „ d /m )/^2^: j at r = T end = 11.6 correspond- 
ing to the onset of the preheating stage y = y en d. 

Short time scale behavior: Fig[6]shows the evolution 
of ipf until t = 50. The horizontal axis used 
here is the dimensionless time r = mot. Clearly, 
the axion masses are all different in the Marcenko- 
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FIG. 6: The evolution of the term J^. ipf that couples to the 
matter field. The time is as in Fig0 Oscillations are less 
sharp than in the equal-mass case, Fig(4] 



Pastur case, resulting in a somewhat different evo- 
lution of the J2i <Pi term from the equal-mass case. 
We can see that the oscillations are more obtuse 
compared to Fig [4] this is a consequence of dephas- 
ing of the axion oscillations owed to relative mass 
differences. As the time-dependent mass term os- 
cillates, some resonance effects for the dynamics of 
Xk are expected. This is indeed the case, at least 
to some extent. Fig(7] (a) shows the time evolu- 
tion of the matter field mode function (g = 10~ 3 
and k/arriQ — 6.0 at r = T em i, as in the previous 
section). The temporal enhancement of the ampli- 
tude (which is clearly seen for small r but becomes 
weaker for large r) is caused by parametric reso- 
nance with (the collective behavior of) the axions. 
In contrast to the equal-mass case, the amplitude 
of Xk decreases on average; the resonance effect is 
not strong enough to resist dilution due to cosmic 
expansion, even when the coupling constant is as 
large as g ~ 10~ 3 . 

One can separate out the effect of cosmic expansion 
by looking at the comoving field Xk = a 3 ^ 2 Xk- The 
equation of motion for Xk is 



-7 + .9 



~(2# ■Ml 2 ) 



X k = 0, (55) 



where the last term in the square bracket is pro- 
portional to the pressure and is very small during 
reheating. FigJT] (b) shows the numerical plot of 
the evolution of Xk- We can see that the peaks 
occur when ^\ ipf of Fig[S] reaches local minima, 
and dJk/^l becomes large (see Fig|8]) , i.e. when 
the system becomes less adiabatic; this is charac- 
teristic of parametric resonance. In contrast to the 
equal-mass case, the minima of the mass term do 
not approach zero due to dephasing, caused by the 
relative mass differences of the axion fields. This 
yields relatively small d>k /w 2 , an< ^ ma kes preheating 
inefficient. 




FIG. 7: (a) The evolution of the mode function of the mat- 
ter field Xk, in JV-flation using the MP mass distribution. 
The coupling is g — 1CP 3 and the wavenumber is chosen 
as k/amo = 6.0 at r = r en d = 11.6. (b) The evolution of 



Xk 



,3/2 



Xk for the same parameters. There are wiggles 



in the oscillation amplitude (these are evident for small r 
and become smaller for larger times) indicating some effect of 
parametric resonance. This resonance is, however, not strong 
enough and the amplitude of Xk decays on average. 



One can see that the amplitude of Xk exhibits 
power-law like growth on average. This growth 
does not necessarily mean production of particles, 
since it is mainly due to redshift [25|. When k 
is large the power Xk ~ a 1 tends to 7 ~ 0.5, 
which is understood as follows: in the mass term of 
(|55[) . k 2 /a 2 is dominant and it stays dominant since 



i- 2 , H ~ t- 2 and H 2 



as 



a ~ t 2 / 3 during preheating. Then ([53]) becomes 
Xk + C't~ 4 / 3 Xk = for some constant C, which is 
exactly soluble in the form of Xk ~ t a F(t) where 
F(t) is a fast oscillating function; discarding the 
decaying solution we find 7 = 3a/2 = 0.5. For 



smaller k, Xk grows faster than 



,0.5. 



for k w 



we find 



,0.75 



numerically. 



Long time scale behavior: in the equal-mass model 
(see section ITVB|) . with a large enough value of g, 
the resonance parameter is q 3> 1 and resonances 
arise for reasonably long time scales, specifically, 
up until r w 4000 for g = 10~ 3 (ignoring backreac- 
tion). Similarly, in the MP case, even though there 
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-0 .03 



FIG. 8: The adiabaticity parameter lu^/cuI, for r < 50. The 
slight negative shift is due to the cosmic expansion. 



is no well-defined g-parameter, resonances can en- 
sue during short time intervals for large r, again 
assuming a similar large coupling g. In this case, 
the collective behaviour of the axions is crucial and 
the adiabaticity parameter iu^/ujj, shows a rather 
complicated behaviour (see Fig[5] (a)). Since the 
mass differences between the neighbouring axions 
in (fTTj) with M = 1500 is typically of order of 
Am 2 ss O(10 -2 ) x m 2 ,, once dephased, the axions' 
collective oscillations return to near-coherence in 
time scales of order At w 0(1O 2 ), causing beats in 
the effective mass for Xk (see Fig|9] (b)). In Fig. 
[TU1 we show the evolution of Xk until r = 2000 
(a) , and the evolution of the comoving occupation 
number calculated for Xk (b). In this example, 
there is some amplification due to parametric res- 
onances around r w 450. On these time scales, for 
g > 10 -3 and small k, we find the occasional am- 
plitude enhancement of a few orders of magnitude. 
For larger wavenumbers (k/amo > 10 4 at r = T en d) 
we find somewhat different behaviour of rife. The 
overall amplitude of uj^ /w 2 becomes smaller but the 
spikes at large t remain. Consequently, the bursts 
at r w 450 disappear and the late time dynamics is 
dominated by a random- walk like behaviour. These 
resonance effects are, however, not frequent or long 
enough to dominate preheating. 

To summarize, we saw that preheating of a single mat- 
ter field is not due to explosive particle production in 
TV-Hation; even though there is some amplification, it is 
too weak in small time scales and not very frequent in 
large time scales, to compete with the dilution due to 
Hubble expansion. We have also studied parameters not 
presented above, including larger values of the coupling 
g; for g — 3 x 10~ 3 the resonance is barely sufficient to 
compete with the Hubble expansion. For the MP param- 
eter (3 = 0.7 and 0.9, we found similar results (inefficient 
resonance). The physical reason for the suppression of 
parametric resonance can be understood as follows: the 
axions are all out of phase, averaging out each other's 



^k/^k ( a j Adiabaticity parameter 
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FIG. 9: (a) The long time scale behaviour of the adiabaticity 
parameter Lbk/^l- 0>) The sum of axions' squared amplitudes 
'^2 i ifi, for the time scale r = 200 to 600. The coupling g and 
the wavenumber k are the same as in FigJTJ 



contribution, so that the driving term oc Y) vf in the 
equation of motion for Xk does not provide a coherent 
oscillatory behavior that is needed for efficient paramet- 
ric resonance. Hence, instead of an exponential increase, 
we observe power-law like behavior Xk ~ a 7-3 / 2 in time 
scales t < 200, where 7 is typically between 0.5 and 0.75 
for the parameter region we have studied. For longer 
time scales r w 200 ~ 2000, we find occasional parti- 
cle production, although these are not strong enough to 
dominate preheating. This conclusion differs from the 
common lore, namely, that parametric resonance effects 
are crucial for preheating [J. 



V. CONCLUSIONS 

In this paper we studied the late time dynamics of 
jV-flation, a string motivated realization of assisted in- 
flation, assuming the Marcenko-Pastur mass distribution 
(arising from random matrix theory) and equal-energy 
initial conditions at the onset of slow roll inflation. We 
provided analytic and numerical calculations of the in- 
termediate phase after the slow roll conditions are vio- 
lated for heavy fields, but before preheating commences. 
We find that the majority of the energy at the onset of 
preheating is carried by the axions with light masses, be- 
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FIG. 10: (a) Long time behavior of Xk> exhibiting short lived, 
weak resonances around r « 450. The choice of parameters is 
the same as in Fig[7] (b) The comoving occupation number 
rife calculated for Xk. 



might change the scenario. For instance, it is argued in 
[33l [3^ that the oscillations of multiple inflatons (with 
irrational mass ratios) can enhance drastically the de- 
cay rate (Cantor preheating). This argument is based 
on two pillars: first, theorems in s pec tral theory indicate 
that stability bands vanish [33|, [34, 38] in the case of more 
fields whose masses are not related by rational numbers. 
Second, numerical evidence in two field models indicate a 
slight enhancement of particle production for well chosen 
parameters 38] . In the latter study of Cantor preheating, 
dephasing of fields is unimportant since only two fields 
are considered. Nevertheless, an enhancement effect like 
in Cantor preheating cannot be excluded in 7V-flation, 
especially if only a few axions couple to a given matter 
field. 

Besides Cantor preheating, we would like to comment 
on yet another effect. It has been shown in [54|, [5f| that 
noise on top of an oscillating driving force can also en- 
hance resonant particle production 3 . This phenomenon 
could also occur in multi-field inflation, if preheating is 
dominated by one or two fields: the oscillations of the 
many other fields would then act similar to noise, poten- 
tially enhancing preheating. 

Further, we have ignored back reaction and rescatter- 
ing during preheating, since modifications due to these 
two effects should be minor. Explosive x~P ar ticle pro- 
duction due to parametric resonance is irrelevant in J\f- 
flation, as argued above, so that (x 2 ) remains small; in 
addition, their inclusion would only diminish resonance 
effects further. 



cause ~ 90% of the energy is carried by only ~ 10% of the 
fields. Thus, only these light fields need to be taken into 
account during preheating. To study preheating, we cou- 
pled a single massless bosonic matter field x to the axions 
<Pi, assuming the same coupling constant g 1 between \ 
and <pi. Within this setup, we solved for the evolution of 
the matter field numerically, including the expansion of 
the universe, and found power-law like behavior in short 
time scales and occasional, not very frequent resonance 
amplifications in long time scales in the parameter region 
that would give rise to stochastic resonance in single field 
models. In particular, the growth of the matter field is 
generically not strong enough to resist redshifting due to 
cosmic expansion. As a result, the old theory of pertur- 
bative preheating (see e.g. [25|]) applies to this scenario 
and not parametric resonance models. The outcome is 
desirable for the model, as there is no danger of produc- 
ing unwanted relics. The prediction of this model is hence 
rather different from the accepted view that parametric 
resonance effects are crucial for preheating Q . 

The analysis presented in this paper is dependent upon 
several assumptions, such as the chosen matter content 
and coupling constants. In particular, we considered only 
one matter field coupled to the whole spectrum of ax- 
ion fields with the same coupling strength g. Although 
we believe this choice to be reasonable and the model 
to be rather generic, dropping some of the assumptions 
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APPENDIX A: FURTHER ANALYTIC 
ESTIMATES 

In section UlI Al we developed a semi-analytic solution 
(J) which underestimates the numerical values of the po- 
tential energy W (see Table [J). Here we present a second 
analytic approximation, which we call solution (II) , that 



3 However, see 
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gives an upper bound for W after the slow roll condition 
for the heaviest field is violated, to check the numerical 
solution in Table Q] In addition, both analytic solutions 
can be used for arbitrarily large numbers of fields that 
would not be tractable via a brute force numerical inte- 
gration. 

The basic idea consists of holding fixed heavy fields 
as soon as the corresponding 77 becomes of order one: 
first, we take the continuum limit so that we can make 
use of the Marcenko-Pastur law for the continuous mass 
variable 1 < x < £. Second, we partition this interval 
into M bins according to a simple rule and denote the 
upper boundaries of bins with Xa, A = 1,...,M, so 
that Xm = Third, whenever t]a (corresponding to 
some Xa) becomes of order one, we hold fixed all fields 
with masses in the A'th bin. Naturally, one recovers the 
full microscopic model if one takes M, = J\f and uses the 
Marcenko-Pastur law as a rule for choosing the bins so 
that Xa = fh 2 A /rriQ. 

Taking M. < J\f leads to a coarse-grained model which 
is more tractable, but one pays the price of having a 
larger W . This approximation is justified as long as the 
energy left in the heavy fields is small compared to the 
energy in the light fields; M. ~ 50 suffices for the range 
of y- values that we are interested in 4 . 

We now proceed to compute Wn , 07/ and Nji as out- 
lined above. We assume first a partition {X\ 1 . . . ,Xm} 
of the interval 1 < x < £ and denote with Ya the values 
of y where = 1 (note that Ya < Yg if A > B). 

If we further denote the energy Wn that is valid in the 
range Ya < y < Ya-i with Wa, we can calculate the 
corresponding Ya as the solution to 



Wa(Ya) = mlX M - A +i 



(Al) 



starting with W\ = W>. Note that Y 1 = yx from (|3"T|) . 
as it should. We can then compute Wa for A > 2 to 



W A (y) 

N 2 *2 



(A2) 



(y x ) 



Xm—a+i 
1 



A-l 

Ere 




Similarly, if we denote with a a the effective field which 
is valid in the range Ya < y < Ya-i (so that <j\ = 07), 
we arrive at 



a A (y) = (ta-i(Ya-i) 

"Pa I ~ 

y 



+ 



ds (A3) 



and finally the number of e-folds becomes (with Yq = 1 

W A (s 



N A (y) = Na-x(Y A -i) + 



-ds . 



(A4) 



In an appropriate large .M-limit the above approxima- 
tion becomes independent of the partition, which is of 
course our aim. We would like to use Wji up to when 
Wi is no longer a viable lower bound for the true energy 
W, that is, until a ~ cri(y e nd)- This is possible by tuning 
X 1 such that a n {Y M ) w (Ji{y en d) 5 - 

That way, the energy ratio of heavy to light fields be- 
comes 



R 



j^W n (Y M )-{Y^) 



Wlight 



Ml 



1 



(A5) 



which has to be smaller than one (see Table UII|) . so that 
we can trust our approximation. 



M 


Xi 


(tii(Ym) 


Wu(Y M ) 


yv n (Y M ) 


NliVend) 


R 


200 


2.90 


5.40 


2.91 


1.47 


1.07 


0.97 


400 


2.15 


8.32 


2.15 


1.24 


1.01 


0.74 


1500 


1.75 


17.57 


1.75 


1.17 


1.00 


0.53 



TABLE III: Semi-analytic solution II for M = 200, 400, 1500, 
to be compared with the slow roll result I from table U Here, 
R is the ratio between the potential energy of the heavy (held- 
fixed) fields to that of the light (dynamical) ones. We use a 
mass splitting into M. — 50 bins. X\ is chosen such that 
vii(Ym) ~ cri(y e nd)- The approximation (//) approaches 
the slow roll result for increasing J\f. The numerical results 
in Table U lie nestled between the two approximations. 



We compare (I) and (II) solutions in Table IIII1 where 
we also vary the number of fields. The solutions approach 
each other in the large Af limit. Henceforth, the numer- 
ical solution is well approximated by either one in the 
case of A/"-flation, where we deal with thousands of fields, 
and consequently, we are justified to use the slow roll 
approximation to set the initial stage for preheating. 



4 Note, generically crj(y) ^ au(y) because firstly, a number of 
fields are artificially held fixed and no longer contribute to the 
path length a, and secondly, the total potential energy is bigger 
so that light fields evolve slightly slower. Only small corrections 
result, since fixed fields are already near the minimum of their 
potential, not contributing much to a anyhow. Moreover, we 
demand R < 1. Consequently, we have yn ~ yj and ajj ~ ctj. 



5 We choose Xi as large as possible so that R < 1, while 
keeping X\ small enough to ensure that the solution (II) re- 
mains applicable up until 07 leaves slow roll; thus we demand 
fir (Xm ) ~ "'/(j/end): leading to X\ = 1.75 for N = 1500. Simul- 
taneously, to distribute the remaining bins, we choose (M — 2)/2 
narrow bins from X\ to X M /2 ~ H (the MP distribution peaks 
in that region), followed by larger bins up until X^a = £ fa 34. 
For Af = 1500 even M < 50 yield results insensitive to the chosen 
partition. 
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